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ABSTRACT 

GX 339-4 has been one of the key sources for unravelling the accretion ejection cou- 
pling in accreting stellar mass black holes. After a long period of quiescence between 
1999 and 2002, GX 339-4 underwent a series of 4 outbursts that have been intensively 
observed by many ground based observatories [radio, infrared(IR), optical] and satel- 
lites (X-rays). Here, we present results of these broad-band observational campaigns, 
focusing on the optical- IR (OIR)/X-ray flux correlations over the four outbursts. We 
found tight OIR/X-ray correlations over four decades with the presence of a break in 
the IR/X-ray correlation in the hard state. This correlation is the same for all four 
outbursts. This can be interpreted in a consistent way by considering a synchrotron 
self-Compton origin of the X-rays in which the break frequency varies between the 
optically thick and thin regime of the jet spectrum. We also highlight the similarities 
and differences between optical/X-ray and IR/X-ray correlations which suggest a jet 
origin of the near-IR emission in the hard state while the optical is more likely dom- 
inated by the blackbody emission of the accretion disc in both hard and soft state. 
However we find a non negligible contribution of 40 per cent of the jet emission in the 
V band during the hard state. We finally concentrate on a soft-to-hard state transition 
during the decay of the 2004 outburst by comparing the radio, IR, optical and hard 
X-rays light curves. It appears that unusual delays between the peak of emission in 
the different energy domains may provide some important constraints on jet formation 
scenario. 

Key words: accretion, accretion discs - binaries: general - ISM: jets and outflows 
- infrared: stars - radio continuum: stars - X-rays: binaries - stars: individual: GX 
339-4. 



1 INTRODUCTION 

Black hole X-ray binaries (BHXBs) are known to be power- 
ful multi wavelength emitters spending most of their life in 
quiescence and undergoing sporadic outbursts during which 
their X-ray luminosities can increase up to a factor 10® com- 
pared to quiescent level. During these outbursts, BHXBs 
display a sequence of three main X-ray states, which are 
defined by their X-ray spectral and timing properties. The 
hard state, is usually observed at the beginning and the end 
of a 'typical' outburst. The X-ray spectra are dominated by 
non-thermal emission in the form of a power-law extending 
up to hard X-rays. The physical origin of this emission is 
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still debated, but the common models imply a corona of hot 
plasma surrounding the compact object or the base of com- 
pact j ets (Zdziarski ct al. 1998; Nowak, Wilms fc Dove 2002] 



Markoff, Nowak fc Wilms|2065| ". The soft state is dominated 



by thermal emission from the accretion disc and displays a 
weak power-law component with photon index softer than 
during the hard state. The transitional states between the 
hard state and the soft state can be classified as various 
instances (hard or soft) of the Intermediate state and has 
spectra with hardnesses in between the hard and the soft 
states. For a complete description of X-ray states properties 



and classification see e.g. Homan fc Belloni| ( 2005 1 , McClin- 



tock fc Remillard^ (2006, ) and |Fender| ( |2006| 



Among the canonical states, it is perhaps the hard state 
that has attracted most attention in recent years where radio 
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observations have highlighted analogies with the flat spec- 



tra of low-luminosity active galactic nuclei (AGN; Falcke & 



Biermann 19961. The radio spectrum from GX 339-4 dur- 



ing its hard state can be explained as self-absorbed syn- 



chrotron emission from compact jets (Corbel et al. 2000 



Fender 2001 1, similar to the case of flat spectrum AGN 
( Blandford fc Konigl|1979 1 . Regular radio observations have 
shown that the compact jets of GX 339-4 are quenched in 
the soft state ( [Fender et al.|1999| [Corbel et al.|2000| ). Simi- 
lar properties have now been found in a growing number of 
BHXBs (the jet has even been resolved in Cyg X-1 



Stir- 



ling et al.|[200Tj and GRS 1915-H0 5, [Dhawan, Mirabel fc 



Rodriguez 2000 Fuchs et al. 2003), thus suggesting that 



compact jets are ubiquitous in BHXBs during the hard state 



(Fender 20011. In addition to being responsible for the ra- 
dio emission, the compact jets may dominate the infrared 
emission and could also contribute significantly in the opti- 
cal domain (Corbel et al. 2001; Jain et al. 2001 Corbel & 



Fender 2002 



Homan et al.,,2005 Russell et al. 20061. The 



compact jets may also produce some of the observed X-ray 
emission ( Markoff et al.|2003 2005 1. Thus, the compact jets 
may dominate large fractions of the spectral energy distribu- 
tion (SEP) of BHXBs ( Render 200^ [Corbel fc Fender|2002 
GaUo, Fender fc Pooley^2003i ) 



A good way to assess the contribution of jets at high en- 
ergy and disentangle its emission components, is to perform 
broad-band observations of BHXB simultaneously at radio, 
optical-infrared (OIR) and X-ray frequencies, and study the 
correlations between these frequency domains. Radio/X-ray 



(Corbel et al. 2000 2003 GaUo et al. 20031 and OIR/X 



ray ( [Homan et al.|2005[ [Russell et al.|2 006) flux correlations 
have already been found for BHXBs in their hard state indi- 
cating a strong coupling between accretion and ejection pro- 
cesses. GX 339-4 was the flrst galactic black hole for which a 
strong non-linear correlation between radio and X-ray emis- 



sion was observed in the hard state (Corbel et al. 2000 



2003|). This correlation has been extended to other galac- 
tic black holes (|Gallo et al.[[2003[) and even AGN ([Merloni 



Heinz fc di Matteo[[200 3': 'Falcke , Kording fc Markoff[[200'^ 
Kording, Falcke fc Corbel 2006). In this work, we present 



the detailed evolution of the IR/X-ray flux correlation of 
GX 339-4 over its last flve years of activity. In section 2, we 
describe the observational data set and our analysis method. 
After a brief description of the recent activity of GX 339-4, 
section 3 presents the results of our broad-band study and 
we discuss and propose possible interpretations in section 4. 
Our conclusions are summarised in section 5. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Optical and Infrared 

Optical and IR photometry of GX 339-4 were conducted 
between UT 2002 January 22 and 2007 October 02 (MJD 
52296-54375). Approximately 840 observations were made 
with the SMARTS (Small and Medium Aperture Research 
Telescope System) which currently uses ANDICAM camera 
on the 1.3m CTIO telescope. ANDICAM takes simultaneous 
optical and IR images over a variety of band passes. In the 
observations reported here, we used Johnson-Kron-Cousins 
V and / fllters ( Bessell, Castelli fc Plez|1998 l and standards 



CIT/CTIO J and H fllters ( |EUas et al.|[T982| . The optical 
and IR light curves for the all period are presented in [Buxton[ 
fc Bailyn (20071. We converted the observed magnitudes m 
into spectral flux densities F,^ using the optical extinction 



Av = 3.7 ± 0.2 (Zdziarski et al. 19981 and the extinction 



law of Cardelli, Clayton fc Mathis (19891. The uncertainty 



on the optical extinction has been propagated in the derived 
errors on the OIR fluxes and dominates over the intrinsic 



2.2 X-ray 

2.2.1 Data reduction 

The X-ray observations were performed with the Propor- 
tional Counter Array (PCA) and the High Energy X-ray 
Timing Experiment (HEXTE) on-board the Rossi X-ray 
Timing Explorer (RXTE). We analysed public data taken 
between 2002 January 29 (MJD 52303) and 2007 October 6 
(MJD 54379), corresponding to a total of 622 pointed obser- 
vations. Spectra were produced for each observation using 

HEASOFT V6.4. 

For the PCA, only data from Proportional Counter Unit 
2 (PCU2) were used for the analysis in this work, as it is the 
only operational unit during all observations and is the best 
calibrated detector out of the 5 PCUs. PCA spectra were 
extracted from the standard 2 mode data and a systematic 
error of 0.6 per cent was added. 

For the HEXTE, starting on 2005 December 12, we only 
used cluster B data as cluster A started to lose its back- 
ground measurements capabilities. HEXTE spectra were 
produced from the standard mode archive data and were 
dead-time corrected. 



2.2.2 Spectral analysis and state classification 

The PCA (3-25 keV) and HEXTE (20-150 keV) spectra 
of each observation were fltted simultaneously in XSPEC 
Vll.3.2 using an overall normalization constant that was 
allowed to float for cross-calibration purposes. To fit the 
spectra, we used several combinations of the following mod- 
els: a (cut-off) power-law (cutoff pi or powerlaw), a multi- 
temperature disc blackbody model (ezdiskbb), a Gaussian 
emission line at 6.4 keV (gaussian), a smeared absorption 
edge (smedge), and an absorption component (wabs). The 
hydrogen colu mn density, Nh, w as fixed to a value of 5 x 10^^ 



atoms cm (Kong et al. 2000). At fainter fiux, when GX 



339-4 was not significantly detected with HEXTE, fits were 
made to the PCA spectrum only. We finally obtained an 
average reduced of 0.96 with a minimum of 0.35 and a 
maximum of 1.87. Unabsorbed fiuxes were measured in the 
3-9, 9-20 and 20-100 keV energy ranges (if HEXTE data 
were used). 

We classified the observations into the various X-ray 
states (hard state, intermediate state, soft state) using the 
dates of the state transitions provided by previous spectral 



and timing studies (Belloni et al. 2005 



Smith, Homan & 



Bezayiff 2005: Bellomet al.|[2006, Kalemci et al.||2007| [Del 
Santo et al.[[2009l. When not available, we determined the 



state using the value of the power-law photon index (hard 
state: F < 2.1, soft state: F > 2.1) and the evolution of the 
hardness ratio. 
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2.2.3 Galactic ridge emission 
At very low flux level (< 10~" 



ergs 



, X-ray emission 



from the Galactic ridge ( Valinia & Marshall 1998 1 can signif- 



icantly contaminate the estimated GX 339-4 fluxes. In order 
to subtract this contribution, we used quasi-simultaneous 
Chandra and RXTE /PCA observations of GX 339-4 taken 
on MJD 52911 during a quiescent state of GX 339-4 ( |Gallo, 
Fender fc Corbel|2003 1, to estimate the Galactic ridge flux, 
Pgr. We analysed the seven PGA observations over this pe- 
riod and combined the spectra to extract an average 3-9 keV 
flux. The high spatial resolution of Chandra provides an es- 
timated real flux for GX 339-4 and is then subtracted from 
the averaged PGA flux. We flnally obtained, in the 3-9 keV 
range. For = (4.2±0.5) x 10"^^ ergs~^cm~^. Consequently, 
we subtracted this value to all 3-9 keV fluxes. 



2.3 Radio 

Since 1996, we have been performing regular radio observa- 
tions with the Australia Telescope Compact Array (ATGA) 



several times a year and during each outburst (e.g. Fender 
et al.|1997| [Corbel et al .'2000, 2003). In this paper, we focus 



on the radio data up to 2005 that is presented in Corbel et al. 
(in preparation). The ATGA synthesis telescope is an east- 
west array consisting of six 22 m antennas with baselines 
ranging from 31 m to 6 km. The continuum observations 
have been mainly performed in two frequency bands (with 
a total bandwidth of 128 MHz for 32 channels), usually at 
4800 MHz and 8640 MHz simultaneously. All ATGA results 
are summarised in Corbel et al. (in preparation) , which also 
includes further information on the radio analysis as well 
as details on the evolution of the radio/X-ray correlation of 
GX 339-4. 



3 RESULTS 

3.1 Five years of activity 

Fig. [l] represents the RXTE/ ASM light curve between 2002 
and 2007 (upper panel) and the light curves in the 3-9 keV 
and 20-100 keV ranges for the 622 analysed RXTE observa- 
tions (lower panel). During this period, GX 339-4 underwent 
a series of four outbursts separated by periods of quiescence. 
The first one and one of the most luminous outburst of those 
recently observed for this source began in 2002 and ended 
in 20 03 ( [Miller et al.|2004[ [Belloni et al.|2005[ [Homan et aT 
2005 ). During the 2002 hard-to-soft state transition, a bright 
radio flare was observed ( [Fender et al.|2002[ ), which later on 
lead to the formation of a large-scale relativistic jet ( Gallo 



et al.|2004[ ). GX 339-4 entered a s econd outburst in 2004 and 
went back to quiescence in 2005 ( [Belloni et a"L][2006l [Miller[ 



et al.[[2006| [Joinet et al.|[2007[ ). During this outburst. Cor 



bel et al. (in preparation) show the presence of two parallels 
tracks in the radio/X-ray correlation between the rising and 
decaying phases of the outburst. During the following 'mi- 
nor' outburst in 2006, GX 339-4 remained in the canonical 
hard state and reached, in the 3-100 keV band, a maximum 
unabsorbed flux of 2.75 x 10'^ ergs"^ cm~^ (See Fig. [l| . 



outburst. Miller et al. (20081 reported the detection of rela- 
tivistically broadened iron emission line in the X-ray spectra 
which would suggest a black hole spin parameter close the 
maximal value. In the following, we will refer to these out- 
bursts by their starting year, namely, 2002, 2004, 2006 and 
2007. 

3.2 OIR/X-ray correlations 



et al.|2001||Gorbel et al.|2001| |Corbel & Fender||2002| |Bux- 


ton & Bailyn|2004 


Homan et al.[[2005 


Russefl et al.|2006l, 



The last outburst occurred in 2007 (Tomsick et al. 2008 



Del Santo et al.|2009| ) and displayed a luminosity similar to hole mass M = 5.8 Mq and a distance to the source D 



provide a basic picture of SED, with the compact jets emis- 
sion dominating from radio up to near-IR (NIR), the accre- 
tion disc dominating the optical-to-UV range and a power 
law component (corona or base of the jets) dominating the 
soft-to-hard X-rays. 

In order to study the connections between these three 
components, we plot our quasi-simultaneous {At 1 day) 
infrared H band (Figs. [2^, [3j [4| and optical V band flux 
density (Fig. ^p) against the 3-9 keV X-ray flux for the four 
outbursts. In Fig. [2^, we show the evolution of the IR/X-ray 
correlation over the four outbursts with the data separated 
by X-ray state (hard state, intermediate state, soft state). 
In Fig. [3j we highlight the difference between the rising and 
decaying phases of each outburst in order to search for ex- 
istence or not of a 'parallel track' effect (see below) in the 
IR/X-ray correlation. In Fig. [4j the same data set is sep- 
arated by outbursts to compare their evolution with time. 
Finally, the optical/X-ray correlation separated by X-ray 
state is represented in Fig. [2j3. 

3.2.1 Infrared vs. X-ray 

Fig. [2|i shows two distinct patterns for the hard state and 
the soft states that are connected by the intermediate-state 
data. In the hard state, IR and X-ray emissions are strongly 
correlated over four decades in X-ray flux. We note a possi- 
ble deviation from a straight line: there seems to be a break 
in the correlation to a steeper slope for data points at X- 
ray fluxes below ~ 5 x 10^^^ ergs~^cm~^. To determine the 
significance of this break, we fit the hard-state data in log 
space with a simple power law and a broken power law. The 
broken power law improves the quality of the fit with a of 
584 (for 130 degrees of freedom (d.o.f.)) compared to a of 
848 (for 132 d.o.f.) for the simple power law. In addition, the 
confidence intervals at 90 per cent of the broken power law 
slopes bi and 62 do not present any overlap (0.64 ^ foi ^ 0.74 
and 0.47 55 &2 ^ 0.49). We performed an F-test to compare 
the two model taking the simple power law model as the null 
hypothesis. We obtained a probability of 3 x 10~^^, which 
indicates that the deviation is indeed signiflcant. We note 
however that both models give a reduced signiflcantly 
greater than 1, so there must be some source of intrinsic 
variability away from the power law relation that dominates 
over the observational errors. We found that the break oc- 
curs at Fbreak = (3.3li;o) >< 10" " erg s" ^cm"^ in the 3-9 
keV range. The corresponding bolometric (3-100 keV) flux 
has an average value of 1.1 x 10"^" ergs~^cm~^. This gives 
a break luminosity Lbreak ~ 10"'^ I/Edd if we assume a black 
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the one reached in 2002. During the initial hard state of this kpc ( Hynes et al.|2063 2004 1 
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52500 53000 53500 54000 

Time (MJD) 



Figure 1. Upper panel: RXTE/ASM light curve of tiie four GX 339-4 outbursts between 2002 and 2007. Ttie red parts of the light curve 
represent the periods covered simultaneously in X-ray and OIR. Lower panel: RXTE/PCA and RXTE/HEXTE X-ray light curves in 
the 3-9 keV and 20—100 keV ranges corresponding to the set of observations analysed in this work. 




1«10"'^ 1x10"^^ 1x10"'^ IxlO"^" 1x10"^ 1x10"^ 1x10"^ 1x10"'^ 1x10"^^ 1x10"^^ l-IO"^" 1x10"^ 1x10"^ 1x10"^ 

X-ray flux 3-9 keV ( erg s"^ cm'^ ) X-ray flux 3-9 keV ( erg s"^ cm'^ ) 



Figure 2. (a) Left panel: Quasi-simultaneous infrared H-band flux density versus 3-9 keV X-ray flux during the four outbursts. Data 
are separated by X-ray states. Red, green and blue points correspond to hard, soft and intermediate state respectively. Red and green 
continuous lines indicate the fit to the hard state and the soft state respectively. Grey zone indicates the 90 per cent confidence interval 
of the X-ray break flux -Fbroak- (b) Right panel: Quasi-simultaneous optical V-band flux density versus 3-9 keV X-ray flux during the four 
outbursts. Data are separated by X-ray states. Red, green and blue data correspond to hard, soft and intermediate state respectively. 
Red and green continuous lines indicate hard state flt and soft state flt respectively. Both plots have been made with the same scales for 
comparison purposes. 
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Figure 3. Quasi-simultaneous infrared H-band fiux density versus 3-9 keV X-ray flux during the four outbursts. Red and green continuous 
lines indicate the broken power-law fit to the hard state and the simple power law fit to the soft state data respectively. Grey zone indicates 
the 90 per cent confidence interval of the X-ray break fiux -Fb^eak • 



Table 1. Results of the non-linear power law fit (i^oiR = '^^x^ ^^'^ broken power law fit (-Fqir = '^^x ^'^^ ^ ^brcak a^nd 
^OIR = "^^x ^ ^break^ ^ ^break) of the OIR/X-ray Correlations with Fqir in mjy, Fx and -Fbrcak in ergs~^cm~^. 

Errors are given at a 90 per cent level of confidence. 



Data set 


X-ray state 


h bi h2 


a 


-^brcak 


H-band / 3-9 keV 


Hard state 
Soft State 


0.68 ±0.05 0.48 ±0.01 

0.34 ±0.01 


(0.7+g:«) X 10« 
(3.5 ± 0.6) X 10^ 


(3.3tJ;«) X 10-" 


y-band / 3-9 keV 


Hard state 
Soft State 


0.44 ±0.01 
0.45 ±0.04 


(3.5 ± 0.3) X 10^ 
(l.rt?;^) X 10^ 





Likewise, in the soft state, a correlation exists be- 
tween IR and X-ray emissions with a power-law slope of 
h = 0.34 ± 0.01. However, the correlation is weaker in the 
soft state (correlation coefficient: 0.84) compared to the hard 
state (correlation coefficient: 0.96 and 0.97 for the low- and 
high- luminosity portions, respectively). This supports the 
existence of a change in the physical processes involved be- 
tween these two states, as we expect the quenching of the 
compact jets and the thermal emission from the accretion 
disc to dominate the soft X-rays. 

In Fig. [4] we note the similarities in the behaviour of 
GX 339-4 over the four outbursts. The source starts at low 
IR and X-ray fluxes along the lower flux power law correla- 
tion and then joins the upper hard-state track. The fluxes 
rise along this track until the hard-to-soft state transition 
is reached. From this point, the IR flux drops dramatically 
while the source reaches the soft state-track. It then evolves 
more randomly along this path, according to the flux varia- 
tion during the soft state. IR and X-ray fluxes decrease until 



the soft-to-hard state transition where the IR flux strongly 
increases, leading back to the hard-state correlation. The 
fluxes then decline following the same path as the rising 
phase. We note however that only five points (from the 2002 
outburst) at fluxes below inbreak belong to a rising phase. 
This means that we cannot rule out the possibility that the 
source does not follow the steep hard-state track during the 
rising state of an outburst. It could be important to keep 
this in mind for further discussions. 

3.2.2 Optical vs. X-ray 

The F-band/X-ray data set (Fig.[2j3) also exhibits two main 
correlation tracks for the hard state and the soft state. How- 
ever, there are several differences with respect to the IR. If 
we fit the hard-state data with a simple or a broken power 
law as previously, the quality of the fit is not improved by 
the broken power-law model. This indicates that either no 
deviation is present in the F-band/X-ray correlation in the 
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Figure 4. Same as Fig.|3]but with data set divided in outburst. Arrows are plotted to indicate the time evolution during the outburst. 



hard state or it is too weak to be statistically significant. 
We note also that the drop in optical luminosity during the 
hard-to-soft state transition is -weaker. These various differ- 
ences may imply that in the optical bands, several other 
emission components, besides the jet emission, contribute 
or even dominate. As already suggested by previous studies 



(e.g. van Paradijs & McClintock 1995 Homan et al. 2005 



Russell et al. 2006 



and references therein), these compo- 
nents may be blackbody emission of the accretion disc or 
reprocessed emission of the X-rays by the outer parts of the 
disc. 

We finally fit the hard- and soft-state data with the 
simple power-law model and we obtained a correlation slope 
of 6 = 0.44 ± 0.04 and 6 = 0.45 ± 0.04 for the hard state 
and the soft state respectively. All fit results for IR/X-ray 
and optical/X-ray correlations are reported in Table [l] 



3.3 Parallel tracks in the hard state 



We focus here on the phenomenon pointed out by |Russell| 
et al.|(|2007|l for the IR/X-ray correlation of XTE J 1550-564 



during its 2000 outburst (phenomenon called 'hysteresis' in 
this article) and for the radio/X-ray correlation of GX 339-4 
over the 1997-2005 period by Corbel at al (in preparation). 
This is the presence of quasi-parallel correlation tracks corre- 
sponding to the hard-state rising phase and the hard-state 



decaying phase of a same outbursl[^ According to Fig. [s] 
and|4] there is no evidence for this particular phenomenon 
in the IR/X-ray correlations of GX 339-4. Indeed, the differ- 
ent hard state phases (rise and decay) of the four outbursts, 
share the same correlation track. Interestingly, we do observe 
two parallel correlation tracks in the radio/X-ray correlation 
during the 2002 and 2004 outburst of GX 339-4. This im- 
plies that, even if there is strong evidences for a jet origin 
of both the IR and radio emission in the hard state, they do 
not seem to be connected to the X-ray in the same way. 

3.4 Radio - IR connection 

In Fig. |5] we present the radio 8.6 GHz, the IR //-band, 
the optical F-band and the hard X-rays (9-200 keV) light 
curves during the beginning of the 2005 decaying hard state 
(for which we have the best radio coverage). A selection of 
radio-to-OIR SEDs with their corresponding data points in 
the IR/X-ray diagram are then presented in Fig.[6|i and[6]3. 

At the beginning of the decaying hard state, we usu- 
ally observe an increase of non-thermal hard X-ray emission 
( Tomsick, Corbel fc Kaaret|2001 1. In most cases, this is also 



^ As discussed in Corbel et al. (in preparation), the two parallel 
tracks in the radio/X-ray correlation of GX 339-4 do not precisely 
correspond to a single track for the rise and a single track for the 



decay, unlike what is seen in XTE J1550-564 ( Russell et al. 2007 1 
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Radio 8.6 GHz 

H-band 

V-band 

X-ray 9-200 keV 




Fig.|6p), the corresponding data point lies in the transition 
zone between the soft-state and the hard-state correlation 
tracks. The second SED, 11 days later (MJD 53493), shows 
a decrease of the radio flux by a factor ~ 3, and an increase 
in IR and optical by a factor ~ 2 and ~ 0.5 respectively. 
We can see that the H-, J- and /- bands fluxes are now 
consistent with an extrapolation of the radio spectrum by 
a power law with a slope a = 0.17 ± 0.01. The connection 
of the radio and IR spectra would indicate an increase of 
the turnover frequency to values greater than the //-band 
frequencies. The shape of the OIR spectrum still suggests 
a significant thermal contribution from the disc emission. 
On Fig. [6]d (green point), the corresponding IR/X-ray data 
point has reached the main hard-state correlation. The last 
SED (MJD 53502) displays features similar to the previous 
one. The IR spectrum is consistent with an extrapolation of 
the radio spectrum with an inverted power law of spectral 
index a = 0.21 ± 0.02 and the optical ( V,I) spectrum still 
suggests the contribution of the disc emission. These three 
SEDs seem to indicate that the radio to infrared spectra get 
more inverted as the fiux (radio and X-ray) decreases. 



Figure 5. Upper panel: Radio 8.6 GHz (mjy), //-band (mjy), 
y-band (mjy) and X-rays 9-200 keV (10"" erg cm-^ s"!) light 
curves at the beginning of the 2005 decaying hard state between 
MJD 53450 and MJD 53550. The radio flux densities have been 
divided by 2 to clarify the plot. Lower panel: The corresponding 
X-ray hardness ratio as a function of time. 



accompanied by an increase in radio and IR emissions, sug- 
gesting the onset of the compact jets once the source has 



settled down back in the hard state 1 


Corbel et al.|2000 Jain 


et al. 2001 Buxton & Bailyn 2004 


Kalemci et al.||2006|. 



Radio, IR and hard X-rays emissions continue to decrease 
while the source goes back to low luminosities. Fig. [5] is con- 
sistent with this description but we note an interesting delay 
between radio and NIR emissions. Indeed, from the first ra- 
dio observation (MJD 53481), radio and hard X-ray fluxes 
are already decreasing while the IR peak flux occurs about 
12 days later (~ MJD 53493). The evolution in optical is 
similar to the IR, with a flux increase of ~7 mJy in V-band 
and ~8 mJy in //-band. 

Fig.[6^ represents three selected radio-to-OIR SEDs cor- 
responding to the dates of the second radio observation (no 
OIR data were available during the first one), the peak of 
the IR fiux and the last radio observation. The first SED 
(MJD 53482) displays a radio spectrum with a power-law 
slope a = 0.13 ± 0.02 consistent with the typical optically 
thick synchrotron emission of a self-absorbed compact jet 
(for a flux density S^ oc We note that the OIR lies 

significantly below the extrapolation of the radio spectrum. 
As the radio emission originates from the optically thick 
portion of the compact jet, the transition frequency to the 
optically thin regime should be located at lower values than 
the //-band frequencies. However, optically thin synchrotron 
emission cannot be dominant in OIR since the OIR spectrum 
increases with frequency. This rather suggests a contribution 
from the outer parts of the accretion disc (disc blackbody. 
X-ray reprocessing). In the IR/X-ray diagram (red point in 



4 DISCUSSION 

4.1 OIR/X-ray correlations in the hard state. 

4.1.1 Are the IR/X-ray correlations consistent with 
synchrotron self-Compton X-rays ? 

During the hard states of each of the four outbursts of GX 
339-4 between 2002 and 2007, we found a strong correlation 
of the form Lir — a L\ between the infrared //-band lumi- 
nosity, LiR, and the 3-9 keV X-ray luminosity, Lx, with a 
slope &2 = 0.48 ±0.01, which breaks into bi = 0.68 ±0.05 at 
X-ray fiuxes below Fbreak = (3.3+i o) x 10"" erg^^cm-^ 



Previous OIR studies of BHXBs in outbursts (|Jain et al. 



2001] [Corbel et al.||2001[ [Corbel fc Fender||2002| [Buxton fc 



Bailyn||2004| Homan et al.||2005 RusseU et al.||2006| 20071 



have presented evidences for a compact jet origin of the NIR 
emission in the hard state. We also know by spectral analy- 
sis that the 3-9 keV X-ray emission is dominated in the hard 
state by a non-thermal component whose origin is usually 
attributed to inverse Compton emission from a hot plasma, 
a corona surrounding the compact object or the base of the 
compact jet. Our results confirm the tight connection be- 
tween these two frequency domains. But thanks to the ex- 
tensive data set collected during a period of 5 years and 
over four outbursts of a given source, it allows us to con- 
strain the correlation index and to reveal some additional 
behaviour that has not been seen in a smaller data set or in 
a sample combining several sources, and hence, introducing 
more scattering due to uncertainties in distance estimates. 
On the other hand, the results and interpretations presented 
here are only valid for one source, and would need to be vali- 
dated for other sources before extrapolating its universality. 

Regarding the correlation indices in the hard state, we 
know that several models have explained the values found 



in previous IR/X-ray or radio/X-ray studies (e.g. Merloni 
et al.|2003||Markoff et al.|2003| [2005} [Yuan fc Cui|2005[ ). So 

the interpretation presented below is not unique but is an 
attempt to explain in a consistent way the results outlined 
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Frequsncy (1 o'" Hz) X-ray flux 3-9 keV ( erg s"' om'^ ) 



Figure 6. (a) Left panel: Radio to OIR SEDs during tlie onset of the compact jet in the 2005 decay, (b) Right panel: The corresponding 
observations in the IR H-hand vs. X-ray 3-9 keV diagram. The continuous lines indicate the fits to the overall data set as already shown 
in Fig.|2]and|3] 



in the previous section and in particular the indices of the 
power-law fit and the break around ~ 10""^ I^Edd- 



As already suggested by Markoff et al. ( 2005 1, the non- 



thermal hard X-ray emission could arise from inverse Comp- 
tonisation of the synchrotron photons produced in the com- 
pact jets, by the plasma producing the synchrotron emission. 
This synchrotron self-Compton (SSC) process, implies a de- 
pendence of the X-ray luminosity, Lx, with the jet power, 
Qjet (Falcke & Biermann 1996 see also discussion in Corbel, 



Koerding & Kaarct 2008 



for V404 Cyg) as: 



Lx ocQ 



11/4 
jet 



(1) 



From the standard conical jet model, we expect the 
monochromatic luminosity at a given frequency L„ to scale 
with the frequency as cc . If the frequency lies in the 
optically thick (flat or inverted) part of the spectrum, we 
usually observe a = 0-0.15 and if it lies in the optically thin 
part, we expect a = —0.5 to —0.7. If we call i/t the break 
frequency between the optically thick and the optically thin 
portion of the spectrum, we get: 



(2) 



As described in |Falcke fc Biermann| |T"996') and lMarkoff et ah] 



1 2003 1 , it follows from simple analytic arguments that oc 



— 17/12 2/3 

Q^gj and vt oc Q^^^ and hence we get the monochromatic 



luminosity at a given frequency: 



ijet 



Using equation ([TJ, we get: 



Lh oc L\ 



(3) 



(4) 



where Lh is the i/-band luminosity. If we consider the typ- 
ical ranges of values of a given above, we obtain b — 0.48- 
0.52 if the _ff-band lies on the optically thick part (a — 
0-0.15) and b = 0.63-0.68 if it lies on the optically thin part 
(a = —0.5 to —0.7). These values are fully consistent with 
the indices we derived for the hard state: b = 0.48 ± 0.01 for 
the high-flux part and 6 = 0.68 ± 0.05 for the low-flux part. 
As the break frequency varies with jet power {ui, oc Q^J^), 
this suggests that the //-band is located on the optically 



thick part of the spectrum for Lx ^ W ^ L^dd and on the 
optically thin part when Lx goes below ~ 10"^LEdd. The 
results are therefore consistent with an SSC origin of the 
X-rays and it is interesting that this interpretation could 
explain in a consistent way the values of the observed in- 
dices and the presence of the break in the correlation (but 
see section |4.1.3| for an alternative interpretation). Fig. [t] 
illustrates the fact that, for a given variation of the X-ray 
luminosity, the corresponding variation in the IR will be 
more important if it lies on the optically thin part of the jet 
spectrum than on the optically thick part, and hence, will 
give a much steeper slope in the correlation. 



4.1.2 Break frequency 

Under the previous assumptions on the origin of the break 
in the correlation, we can calculate an expression giving 
a rough estimate of the break frequency as a function of 
the monochromatic luminosity L^. As presented above, in a 
standard jet model the luminosity depends on the jet power 

17/12 2/3 

as L[^^ oc 

QjJt and the break frequency as Vb oc Qj^t ■ 



Hence, we find L^ 



17/8 



and we can write: 



17/8 



(5) 



where a is the normalisation factor. The broken power-law 
fit gives a value of the i/-band flux where the break in the 
correlation occurs of Fh — 5.2 ± 1.5 mjy. At this flux den- 
sity, we therefore have i^t ~ vh, with uh = 1-18 x lO^'* Hz 
the effective frequency of the //-band filter. Consequently, by 
converting the //-band flux density Fh into the monochro- 
matic luminosity Lh at the frequency uh, we obtain the 
following normalisation factor for an assumed distance to 
the source of 8 kpc: 



(5± 1.5) X 10' 



D 

8 kpc 



erg s 



-^Hz- 



(6) 



Using equation |[2| and ([5| we can finally express the break 
frequency as a function of a measured monochromatic lumi- 
nosity in the jet spectrum, its corresponding frequency 
f and the spectral index a where this frequency is located: 
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(7) 



with the normaUzation constants stated above. Note that 
the above expression has been derived without any assump- 
tion on the physical origin of the X-ray emission in the hard 
state. Additionally, we would like to point out the work of 



Nowak et al. ( |2005J on the scaling of the break frequency 
with the soft X-ray flux in GX 339-4. The authors infer 
the location of the break frequency by fitting the radio and 
X-ray data by a doubly broken power law which implicitly 
assumes that the soft X-ray emission (3-9 keV) is dominated 
by direct synchrotron emission from the jet. However, since 
we do not adopt the same assumptions, we cannot compare 
consistently their work with our results. 

4.1.3 Is the IR/ X-ray correlation consistent with an 
accretion flow ongin of the X-rays ? 

To explain the fundamental plane of black hole activity, 



Merloni et al. ( 2003 1 tested several classes of accretion flow 
model that could be responsible for the X-ray power law 
emission. Based on the initial work by |Heinz fc Sunyaev] 
(20031, they provided analytical expressions to determine 



the correlation coefficient between radio and X-ray fiux that 
we can use to test these models with our data. We have 
selected the two classes of models the authors identified as 
consistent with their results: the ADAF case and a more 
general class of radiatively inefficient, mechanically cooled 
accretions flows that are modifled by convection or outflows 
and where we expect the X-ray spectrum to be dominated 
by bremsstrahlung emission. 

Note that in the following calculations, we used a power 
law index of the electron distribution p — 2-2.4. For the 
ADAF case with an optically thick synchrotron IR spectral 
indej{^(a = 0-0.15), we obtain b = 0.56-0.62. In the opti- 
cally thin case (a = -0.5 to -0.7) we get b =0.76-0.80. 
The bremsstrahlung dominated case gives slightly higher 
values of b : for a = 0-0.15 we obtain b = 0.66-0.71 and 
for a = -0.5 to -0.7 we get b = 0.88-0.92. 

Consequently, in the frame of an accretion flow origin 
of the X-rays, an alternative interpretation of our results 
would be to consider that the //-band is constantly located 
in the optically thick part of the jet spectrum (a = 0-0.15) 
and that the X-ray emission originates from an ADAF at 
high luminosities that becomes bremsstrahlung dominated 
at low luminosities. This would give 6 = 0.56-0.62 at high 
luminosities, which is roughly consistent with our derived 
index b = 0.48 ± 0.01, and would become b = 0.66-0.71 at 
low luminosities, fully consistent with our b = 0.68 ± 0.05. 

4.1.4 Optical/X-ray correlation in the hard state 



As reported in Section |3.2.2| in the hard state the optical 
emission shows a correlation with the X-rays in the form 
of a power law with slope b = 0.44 ± 0.01. Unlike the IR 
versus X-ray correlation, it does not show any significant 
break in the correlation. Following equations ([5| and (|6| 



^ Note that the spectral index a used by [Merloni et ah] l |2003l l , 
is defined with the following convention: F^, oc which is the 

opposite of the one adopted in this paper. 
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Figure 7. A schematic jet spectrum which illustrates that for a 
given variation of the jet power, the corresponding variation of the 
luminosity will be stronger in the optically thin regime compared 
to the optically thick portion, due to the variation of the break 
frequency. 



and assuming that the synchrotron emission from the jet 
dominates the optical domain, we would expect the break 
to occur at a V-band flux around 130 mjy, which is above 
the maximum flux reached in our data set (Fig. [2j3). Con- 
sequently, if the optical emission is related to the compact 
jets optically thin synchrotron, we would indeed expect a 
simple power law to flt the entire optical/X-ray hard-state 
data with a correlation index b close to 0.68 as seen in IR 
below the break. However the shape of the OIR SEDs and 
the weaker gap in F-band luminosity compared to the IR 
luminosity between the hard and the soft state, rather sug- 
gest that another emitting component than the compact jet 
dominates the optical domain as already suggested by e.g. 



van Paradijs fc McClintock| ( |1995[ ); |Homan et al.| ( |2005 l 
Russell et al.| (| 2006[ ). This component could be the accre- 



tion disc or the companion star and in both cases it can 
be direct blackbody emission or reprocessing of the X-rays. 
[Homan et al.| ( |2005| ) already showed that we can basically 
rule out the secondary star as the dominant source of OIR 
emission in GX 339-4 during outburst. The accretion disc is 
therefore the most likely dominant source in optical in both 
hard and soft state. Consequently, the observed correlation 
between optical and X-ray (3-9 keV) point out the connec- 
tion between the outer parts of the accretion disc and the 
non thermal emission in X-ray which dominates the 3-9 keV 
range in the hard state. Regarding the origin, irradiated or 
not, of the outer disc emission, we can see on Fig.[5]that the 
delay between the hard X-rays and OIR peaks is around 12 
days. If the hard X-rays are indeed the source of irradiation 
and the OIR the reprocessed emission, a 12 days delay is 
highly inconsistent with a reprocessing mechanism. 

Under the above assumptions, the loss of optical emis- 
sion during the hard-to-soft state transition gives an esti- 
mate of the synchrotron contribution in optical of 40 per 
cent in average. 
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Table 2. Expected power law indices, b, in the relation oc L^, 
depending on the emission regime and the disc model. 



Emission regime 


Viscous disc 


X-ray heated disc 


Rayleigh-Jeans tail 


b = 0.25 


b = 0.25 


Flat part 


b = 0.66 


b=l 



4.2 OIR/X-ray correlations in the soft state 

In the soft state, the drop in IR flux is indicative of the 
quenching of the jet emission. The OIR SEDs, that can be 
well fitted with simple power laws suggest that the remain- 
ing IR emission in the soft state originates from the outer 
parts of the accretion disc. In that case, the emission may 
be dominated by thermal emission from a viscously heated 
disc or an X-ray heated disc. 

Let us consider a simple model of a steady thin disc with 
the temperature scaling with the radius as T oc r~", where 
n = 3/4 for a standard disc and n = 1/2 for an irradiation 
dominated disc ( [Frank, King fc R ainc 2002}. Depending on 
the emission regime of the disc spectrum, we have the fol- 
lowing relation between the monochromatic luminosity at a 
given frequency u and the temperature: 



(i) In the Rayleigh-Jeans (RJ) tail: 



(ii) in the 'flat' parlj^of the disc spectrum: 

— 3— — 

If we assume, that Lx oc T^^ oc T*, where Lx is the X- 
ray luminosity of the disc and Tin its inner temperature, we 
obtain: 



L^j^, oc L-^ 



1/4 



and 



-^"tlat °^ ^X 



l/2n 



(8) 



Therefore, we expect L,j^-^^^ oc for a standard disc and 
Lv„ t oc Lx for an irradiated disc. The RJ tail relation 



remains the same in both cases, L„ 



oc lV*. Table ^ 



summarises the expected power law indices depending on 
the emission regime and the disc model. Our derived index 
b — 0.34 ±0.01 for the IR/X-ray correlation in the soft state 
is thus inconsistent with the if-band lying on the 'flat' part 
of a standard or irradiated disc and is only marginally con- 
sistent with the expected relation for the RJ tail. 

However, this index could be explained by considering 
a frequency lying at the transition between the RJ part and 
the flat part of the disc spectrum. Indeed, the correlation 
index 6, expected from the analytical relations given above, 
is inversely proportional to the spectral index a of the disc 
spectrum. Given that this spectral index decays from the RJ 
tail to the 'flat' part, we would expect the corresponding cor- 
relation index to increase from b = 0.25 to & = 0.66 or fe = 1 
for a standard or irradiated disc respectively. Consequently, 
our derived index b = 0.34±0.01 could be interpreted as the 
//-band lying in the transition zone from the RJ tail to the 
'flat' part of the disc spectrum. 



^ Portion of the disc spectrum between the Rayleigh-Jeans tail 
and the Wien cutoff, usually associated with a spectral index 
a = 1/3 in the case of a viscously heated disc. 



This interpretation also implies that the correlation in- 
dex between the IZ-band luminosity and the X-ray luminos- 
ity in the soft state has to be greater than 0.34 since the 
F-band should be further in the transition zone than the H- 
band. This is, indeed, what we observe: b = 0.45 ± 0.04 for 
the optical/X-ray correlation in the soft state. The idea that 
the H- to T/-bands are located in the transition from the RJ 
to the 'flat' part, is also supported by the spectral indices 
derived from the OIR SEDs. Fig. |8] shows the evolution of 
the //-band to F-band spectral index with the optical flux. 
We can see that the spectral index increases with the optical 
flux from ~ 1 to ~ 2. Since we expect the frequency of the 
break between the RJ and the 'flat' part to increase with 
mass accretion rate, and thus, with optical flux, this suggest 
that the OIR bands are 'moving' from the transition zone to 
the RJ tail with the spectral index reaching, at the highest 
fluxes, the expected value a — 2 for the RJ tail. 

However, in the theoretical frame described above, our 
derived indices do not allow us to discriminate between a 
viscously heated or an X-ray heated disc dominating the 
OIR emission in the soft state. This conclusion slightly dif- 



fers from those obtained by Russell et al. (2006, 2008) where 



evidences are provided for an irradiated disc dominating the 
OIR emission in the soft state. These differences arise from 
the way we derive the expected correlation indices between 
OIR and X-ray luminosity in the case of an irradiated disc. 



(2006 



Based on van Paradijs & McClintock (19941, Russell et al 



20081) expect Lqpt oc LJ^ 



for an irradiated disc, 
whereas we expect /^opt oc L'i^'"'~^ using the relations given 
in e.g. Frank et al. (2002 I. This implies that, unlike Russell 
et al. (20061, we cannot conclude on the reprocessed na- 
ture of the optical emission in the soft state. We think that 
further investigations, especially with broad band SED fit- 
ting and time lag studies, are necessary to clearly address 



this issue (see the discussion in e.g. Gierliriski, Done & Page 
|2009[ and references therein). We note, however, that apart 
from this irradiated disc interpretation, our specific results 
on GX 339-4 are consistent with the global interpretations 
of IRussefl etal] ( [20061 [20081 - 

4.3 Radio - IR connection 

Fig. [5] and [6l shows the evolution of the OIR, radio and hard 
X-ray emissions during the formation of the compact jets. 
We first note the ~ 12 days lag of the OIR peak against the 
hard X-rays. This lag has been already observed in several 
BHXBs during outburst decay ( [Kalemci et al.[[2005[ [2006 1 
even if it is usually of the order of ~ 3 days. These has been 
interpreted as the time needed for an optically thin envi- 
ronment (corona, jet base) to be formed before launching a 
jet. We note also that the OIR peak follows the radio by at 
least 10 days. If the increase of the OIR emission is indeed 
due to the formation of the compact jets, this delay implies 
that the onset of the jets synchrotron emission starts at low 
frequency. This is natural since we must have a phase tran- 
sition where the energy density (magnetic, kinetic) of the 
jet has to increase before it reaches its 'stable' structure and 
this would first produce radio emission before OIR as seen 
in Fig. [5] This idea is also supported by the SEDs shown 
in Fig. [6b which suggest an increase of the break frequency 
of the jet spectrum between MJD 53482 and MJD 53493. 
However, an important question to address is whether the 
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V-band flux density (mjy) 

Figure 8. OIR spectral index against V-band flux (mJy). The 
indices are derived from the H- to V-band SEDs during the soft 
state phases only. Dashed line indicate the Rayleigh-Jeans tail 
value of the spectral index. Error bars on the indices and the V- 
band flux are not plotted for clarity purpose and correspond to an 
average value of ±0.2 for the spectral index and ±20 per cent for 
the V-band flux. Note however, that these errors are dominated 
by the uncertainty on the optical extinction and will therefore 
affect all the data in the same way without changing the global 
trend. 



observed delay between radio and OIR is compatible with 
the time scales involved in a jet formation. But this requires 
a detailed model for the creation of a jet which is beyond 
the scope of this paper. 



5 CONCLUSIONS 

Using an extensive broadband data set, we compared the 
connections between X-ray, OIR and radio properties in the 
BHXB GX 339-4. We can summarise our main conclusions 
as follows: 

(i) We first see a strong power law correlation between 
IR and X-ray fluxes in the hard state with the presence of 
a break in the correlation slope below an X-ray luminosity 
around ~ 10~^LEdd- This break and the values of the cor- 
relation indices can be interpreted in a consistent way if we 
consider an SSC origin of the X-rays in the hard state and 
the variation of the break frequency of the jet spectrum. 
However, we note that our results are not inconsistent with 
an accretion flow origin of the X-ray if we consider that 
the X-ray emission originates from an ADAF at high lu- 
minosities that becomes bremsstrahlung dominated at low 
luminosities. 

(ii) Optical and X-ray emissions in the hard state display 
a correlation as well but suggest that the outer parts of 
the accretion disc dominate the optical emission in the hard 
state. 

(iii) In the IR/X-ray correlation in the hard state, we do 
not detect any parallel tracks similar to what is seen for 
XTE J 1550-564 ( [Russell et al]|2006| ) or in the radio/X-ray 
correlation of GX 339-4 (Corbel et al. in preparation). 

(iv) In the soft state, the correlations between the OIR 
and the X-rays indicates also a disc origin of the OIR emis- 
sion and suggest that the //-band and F-band are located 
in the transition zone between the Rayleigh-Jeans tail and 
the 'flat' part of the disc spectrum. 

(v) We compared hard X-ray, optical, infrared and radio 



light curves during a selected transition to the hard state 
and found that the changes in the OIR emission follow the 
changes in hard X-ray and radio. These delays rule out a 
reprocessing mechanism as the origin of the OIR emission 
in the hard state if the hard X-rays are the source of irradia- 
tion, and suggest that the onset of the compact jets emission 
starts by the low frequencies. 

We flnally think that further investigations with good 
enough multi-wavelenght data to follow the evolution of the 
break frequency together with broadband SED fitting will 
allow us to disentangle the dominant emission processes in 
the various energy bands. 
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